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The Dense Plasma Torus Around the Nucleus of an Active Galaxy 

NGC 1052 
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Abstract 

A subparsec-scale dense plasma torus around an active galactic nucleus (AGN) is unveiled. We report on 
very-long-baseline interferometry (VLBI) observations at 2.3, 8.4, and 15.4 GHz towards the active galaxy 
NGC 1052. The convex spectra of the double-sided jets and the nucleus imply that synchrotron emission 
is obscured through free-free absorption (FFA) by the foreground cold dense plasma. A trichromatic 
image was produced to illustrate the distribution of the FFA opacity. We found a central condensation 
of the plasma which covers about 0.1 pc and 0.7 pc of the approaching and receding jets, respectively A 
simple explanation for the asymmetric distribution is the existence of a thick plasma torus perpendicular 
to the jets. We also found an ambient FFA absorber, whose density profile can be ascribed to a spherical 
distribution of the isothermal King model. The coexistence of torus-like and spherical distributions of the 
plasma suggests a transition from radial accretion to rotational accretion around the nucleus. 

Key words: galaxies: active — galaxies: individual (NGC 1052, 0238—084) — galaxies: nuclei 
radio continuum: galaxies — techniques: interferometric 



1. Introduction 

The giant elliptical galaxy NGC 1052, which is also 
classified as a low-ionization nuclear emission-line re- 
gion (LINER) galaxy, holds an AGN in the center (§00 
[cite]cite.Hol997Ho et al.(1997)). If we assume H = 75 
km s- 1 Mpc" 1 and q = 0.5, the redshift, z = 0.0049 (§00 
[cite]cite.Knappl978Knapp et al.(1978)), corresponds to 
a distance of 20 Mpc, and 1 milliarcsec (mas) corresponds 
to 0.095 pc. At radio wavelengths, double-sided jets elon- 
gated by several pc in P.A.~ 65° are seen with an emis- 
sion gap in between (O00 [cite]cite.Claussenl998Claussen 
et al.(1998); §00 [cite]cite.Kellermannl998Kellermann et 
al.(1998)), though it is unclear where the nucleus is lo- 
cated. Water masers are found along the jet, which 
are considered to be excited by shocks into circum- 
nuclear molecular clouds, or amplification of the radio 
continuum emission of the jet by foreground molecular 
clouds (§00 [cite]cite.Claussenl998Claussen et al.(1998)). 
This object shows a convex radio spectrum peaked at 
10 GHz, to be classified as a GHz-Pcakcd Spectrum 
(GPS) source (§00 [cite]cite.ODeal9980'Dea(1998); §00 
[cite] cite. deVriesl997de Vries et al.(1997)). A steep spec- 
trum with cvo = — 1.33 (Si, oc v +a °) at frequencies above the 
peak (§00 [cite]cite.deVriesl997de Vries et al.(1997)) and 
a high brightness temperature of ~ 10 9 K (see subsection 
3.1) indicate synchrotron radiation. 

Two possible hypotheses are proposed for the ris- 
ing spectrum below the peak frequency in GPS 
sources: synchrotron self-absorption (SSA) or free- 
free absorption (FFA). Both can generate convex spec- 



v a ° exp(— Tf v~ 21 ) for FFA, where v is the frequency in 
GHz, t s and Tf are the SSA and FFA coefficients, re- 
spectively, and cvo is the intrinsic spectral index. In 
general, the spectral shape will be a product of both. 
FFA arises from an external absorber in front of the syn- 
chrotron emitter, while SSA is an internal process of the 
synchrotron emission via an interaction between the rel- 
ativistic plasma and the magnetic field. Because these 
spectral shapes resemble each other, precise spectral mea- 
surements are necessary to detect any small difference. 
FFA can generate a steeper cutoff than the maximum at- 
tainable spectral index, a, of 2.5 by SSA. Thus, a com- 
ponent with a > 2.5 could be evidence of FFA. Another 
distinction is related to the difference of the peak fre- 
quencies in twin jets. In the case of SSA, the peak fre- 
quency in the receding jet should be lower than that in 
the approaching jet due to the Doppler shift, if we as- 
sume an intrinsic symmetry On the contrary, a longer 
path length towards the receding jet causes a deeper FFA 
opacity, and thus a higher peak frequency. A sufficiently 
high spatial resolution is crucial to avoid spectral blend- 
ing, as demonstrated by ^|e]cite.Kameno2000Kamcno et 
al.(2000) ([cite]cite.Kameno2000Kameno et al.(2000)). 

Because the FFA coefficient is related to the electron 
density, n c (cm -3 ), and the electron temperature, T c (K), 

by 
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where / LOS dL is the integration via the line of sight 
in pc, it can be a probe for cold dense plasma (§00 
[cite]cite.Pacholczykl970Pacholczyk(1970)). 
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There exists increasing evidence for FFA towards 
active galaxies. One remarkable milestone is the 
discovery of a counter-jet in 3C 84 at the cen- 
ter of the nearby Seyfert/cD galaxy NGC 1275 
(@00 [cite]cite.Vermeulenl994Vermeulen et al.(1994); ©00 
[cite]cite.Walkerl994Walker et al.(1994)). The counter- 
jet showed a strongly inverted spectrum, unlike the ap- 
proaching jet with a steep spectrum, whose feature can 
be explained by FFA due to ambient cold dense plasma 
(see the Appendix). This idea was confirmed by fol- 
lowing VLBA observations by ^D]cite.Walkcr2000Walker 
et al.(2000) ([cite]cite.Walker2000Walker et al.(2000)) 
and VSOP observations by []e]cite.Asada2000Asada et 
al.(2000) ([cite]cite.Asada2000Asada et al.(2000)). 

FFA towards GPS sources was pre- 
dicted by §e]cite.Bicknelll997Bicknell et 
al.(1997) ([cite]cite.Bicknelll997Bicknell et al.(1997)). 
They suggested that photoionization of the interstellar 
medium (ISM) by radiative shocks produces ionized 
gas, which has a significant emission measure and a 
correspondingly high FFA opacity. The discovery of 
FFA towards GPS sources was made by a space- VLBI 
observation towards a GPS galaxy, OQ 208, which showed 
convex spectra towards two radio lobes separated by 10 
pc. The difference in the path length along the line of 
sight between two lobes causes an asymmetry in the FFA 
opacity, with a larger opacity towards the far-side lobe. 
The FFA opacity implied that cold (T e ~ 10 5 K) dense 
(n e ~ 10 5 cm -3 ) plasma associates with the central tens of 
pc (§00 [cite]cite.Kameno2000Kameno et al.(2000)). This 
result showed that multi-frequency VLBI observations 
are crucial to clarify the FFA, which can be a probe 
of matter in the vicinity of AGN. The discovery of the 
FFA absorber was followed by observations towards the 



[cite]cite.Joncs2000Jones 
|00 [cite]cite.Marr2000Marr 



radio galaxies NGC 4261 ( 
et al.(2000)) and 0108+388 
et al.(2000)). 

We emphasize two advantages for the FFA as a probe 
of accreting matter, when the absorber covers the back- 
ground bright core and jets. First, since the absorber must 
be at the near side of the jet, it allows us to figure a 3- 
dimcnsional geometry of the system. Second, the opacity 
provides physical properties of the absorber, such as the 
electron density, n e , and the temperature, T c . NGC 1052, 
the nearest GPS source showing the double-sided jets, is 
one of the best target sources. 

In this paper, we report on simultaneous trichromatic 
Very Long Baseline Array (VLBA) observations towards 
NGC 1052 at 2.3, 8.4, and 15.4 GHz, which produce an 
FFA opacity map. Our aim is to investigate the geometry 
and physical properties of the FFA absorber. 

2. Observations and Data Reduction 

Using the VLBA of the National Radio Astronomy 
Observatory (NRAO), we observed NGC 1052 at three 
frequencies (2.3, 8.4, and 15.4 GHz) simultaneously on 
1998 December 15. The subreflector switched between 
the dual-frequency 2.3/8.4 GHz reflector system and 



the 15.4 GHz feed horn with a 22-min cycle. We 
used 4 channels of 8 MHz bandwidth at 15.4 GHz, 
and allocated 2 channels at both 2.3 and 8.4 GHz. 
The correlation process was accomplished by the VLBA 
correlator. We applied fringe-fitting, data flagging, 
and a priori amplitude calibration in the NRAO AIPS 
(Astronomical Image Processing System). Imaging and 
self-calibration processes were carried out using Difmap 
(0QO [cite]cite.Shepherdl997Shepherd(1997)). We took 
uniform weighting with scaling by errors raised to the 
power — 1. The synthesized beam sizes (FWHM) were 
6.06 x 2.45 mas in P.A.= 3°.92, 2.01 x 0.78 mas in P.A.= 
-1°.91, and 1.03 x 0.40 mas in P.A.= -2°. 53, for 2.3, 8.4, 
and 15.4 GHz, respectively. 

While relative gain errors among the antennas are cor- 
rected though amplitude self-calibration processes, fur- 
ther flux calibration is necessary to obtain certain spectra 
across the observing frequencies. For the purpose of the 
absolute flux calibration, we also imaged four calibrators: 
the BL Lacertae, DA 193, 3C 279, and OT 081. Based 
on a comparison between the total flux measurements by 
the University of Michigan Radio Astronomy Observatory 
(UMRAO) and the NRAO Green Bank interferometer, 
and a summation of the CLEANed flux densities, we ap- 
plied flux scaling for the final results. At 15.4 GHz, for 
instance, the total CLEANed flux densities before an abso- 
lute correction were 3.276, 5.012, 25.769, and 4.172 Jy for 
BL Lac, DA 193, 3C 279, and OT 081, respectively. The 
UMRAO database provides the total flux densities at 14.5 
GHz, 3.458±0.015, 4.763±0.035, 27.69±0.18, and 4.380± 
0.030 Jy for each, as averaged over 2 months centered at 
our observation. We then calculated a correction factor 
(= total CLEANed flux density/UMRAO flux density) of 
0.959 ±0.018, so that we scaled our image by the factor 
and obtained an accuracy of the flux scale of 1.8%. The 
accuracy was derived from a standard deviation of the 
correction factors by the four calibrators. In the same 
manner, we scaled images at 2.3 and 8.4 GHz by correc- 
tion factors of 0.84 and 0.79, respectively. The estimated 
amplitude accuracies are 4.0, 7.2, and 1.8% at 2.3, 8.4, and 
15.4 GHz, respectively. The images are shown in figure [l]. 

3. Discussion 

3.1. Component Identification and Spectrum 

For the purpose of registration, we picked up three dis- 
tinct components (A, B, and C ;see figure [l]). Using the 
task 'IMFIT' in AIPS, we measured the relative positions, 
sizes and flux densities of these components, as listed in 
table 1. In advance of the IMFIT process, we restored 
the CLEAN images with a Gaussian beam of 3 x 1.5 mas 
in P.A.= 0° to match the resolutions. Double-sided jet 
components, A and C, are identified with similar separa- 
tions of 9.5, 7.7, and 6.8 mas at 2.3, 8.4, and 15.4 GHz, 
respectively. An unresolved component, B, is clearly seen 
with a significant flux density of 489 ±9 mJy at 15.4 GHz, 
but becomes fainter 77 ± 6 mJy at 8.4 GHz and too faint 
(< 5.6 mJy) to be identified at 2.3 GHz. 

These components allow us to register three images to 



Plasma Torus around NGC 1052's Nucleus 3 




20 10 -10 -20 

Relative right ascension [mas] 



Fig. 1. Images of NGC 1052 at 2.3, 8.4, and 15.4 GHz, taken using VLBA on 1998 December 15. The synthesized beam sizes 
(FWHM) are 6.06 X 2.45 mas in P.A.= 3°.92, 2.01 X 0.78 mas in P.A.= -1°.91, and 1.03 X 0.40 mas in P.A.= -2°. 53, respectively, as 
shown in the left of each image. The contours start at the ±3<x level, increasing by a factor of 2, where a = 1.246, 2.277, and 0.823 
mjy beam -1 for 2.3, 8.4, and 15.4 GHz, respectively. 
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Fig. 2. Trichromatic color image of NGC 1052 synthesized after restoring the images at 2.3, 8.4, and 15.4 GHz with a 3 X 1.5 mas 
beam, and register by referencing components A and C. The red, green, and blue color channels are allocated for the images at 2.3, 
8.4, and 15.4 GHz, respectively. The origin of all images are referred by the position of component B, which is considered to be the 
nucleus. The blue part covering the nucleus and ~ 5 mas to the west indicates a steeply rising spectrum with a > due to large 
FFA opacities. The red color at both ends of the jets shows a < 0, suggesting optically-thin synchrotron emission. 



Table 1. Parameters of the Gaussian components. 



Comp. 


Frequency 
(GHz) 


Flux Density 
(mJy) 


Relative R. A.* 
(mas) 


Relative Dec.* 
(mas) 


(mas) 


(mas) 


A 


2.274 


969 ± 39 








4.2 


1.2 




8.424 


1800 ±130 








3.4 


0.9 




15.364 


1310 ±24 








3.4 


0.9 


B 


2.274 


<5.6 












8.424 


77.0 ±5.7 


-4.59 ±0.2 


-2.03 ±0.3 


< 1.3 


< 0.7 




15.364 


488.7 ±8.8 


-3.59 ±0.2 


-1.36±0.3 


<0.5 


< 0.4 


C 


2.274 


12.6 ±2.2 


-8.85 ±0.7 


-3.44±1.8 


2.1 


1.4 




8.424 


247.2 ±18.0 


-7.16 ±0.2 


-2.91±0.3 


1.6 


0.6 




15.364 


189. 3± 3.8 


-6.29 ±0.2 


-2.53±0.3 


1.6 


0.6 



*Rclative position with respect to component A. 

tFull width at half maximum (FWHM) of the major and minor axes of the Gaussian component. 



minimize the total difference of their positions. At any 
given frequency, we measured the positions (££,??£) of 
these three components relative to a tentative origin of 
each image, where k and v stand for an index of the com- 
ponents (A, B, and C) and a frequency, respectively. We 
then derived the relative offsets (8£, v ,8rf) between the im- 
ages to minimize the positional residuals, xti defined as 



V 



* 2 (&) + * 2 (C) ^faJO+^W 



, (2) 



where <r(££) and cr(?7fc) arc standard positional errors of 
component k at frequency v. The reference frequency, i/q, 
was put at 8.4 GHz. The r.m.s. errors of the registration 



in (£,!/) were (0.25,0.09) mas between 15.4 and 8.4 GHz, 
and (0.20,0.25) mas between 8.4 and 2.3 GHz, respec- 
tively. The response to the registration errors is discussed 
in subsection 3.4. 

Using image registration, we synthesized a trichromatic 
image (figure |^) allocating red, green, blue channels for 
2.3, 8.4, and 15.4 GHz, respectively, using the beam- 
matched images. The blue color covering component B 
and western ~ 5 mas indicates a steeply rising spectrum 
(a > 0), while the red color at both ends of the elongation 
shows a < 0. 

Figure |^ shows the spectra of the components. Both 
A and C show convex spectra, which can be fitted by 
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Fig. 3. Spectra of components A, B, and C. The flux den- 
sities were measured by a Gaussian model fit using the task 
'IMFIT' in AIPS, as listed in table 1. The errors includes 
both formal errors and amplitude calibration errors. The solid 
lines for components A and C show the best-fit FFA model, 
Sv = Sav ao exp(— Tf i/~ 2,1 ), where S v is the observed flux den- 
sity in Jy, So is the unabsorbed flux density at 1 GHz, ao is 
the intrinsic spectral index, Tf is the FFA opacity coefficient, 
and v is the frequency in GHz. 

optically-thin synchrotron emission obscured by the cold 
plasma via thermal FFA. The spectral index, ao = —1.33 
(©00 [cite]cite.deVriesl997de Vries et al.(1997)), was too 
steep to fit, instead, we obtained ao = —0.65 from the 
best fit. A and C mark the spectral peaks at 4.9 and 
7.6 GHz, respectively, while B must have a peak above 
15 GHz. Although the spectra of A and C seem to be 
similar, the difference of the peak frequencies is signifi- 
cant. The colors of these components in the trichromatic 
image (figure ||) emphasize the difference. The spectral 
indices between 2.3 and 8.4 GHz are a|'3 = 0.5 ±0.1 and 
2.3 ±0.1 for A and C, and > 2.0 for B. Component 
B shows a large spectral index of ag 5 ^ 4 = 3.1 ±0.1 be- 
tween 8.4 and 15.4 GHz. It exceeds the theoretical 
limit for SSA of 2.5 with a uniform magnetic field (©00 
[cite]cite.Pacholczykl970Pacholczyk(1970)). A and C are 
resolved at 15.4 GHz with brightness temperatures of 
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Fig. 4. Distribution of the FFA opacity and a schematic 
model of a torus surrounding NGC 1052. (a): The distribu- 
tion of the FFA opacity, Tf , is drawn in grey scale. The opac- 
ity was calculated by the least-squares fit for the FFA model, 
Sv oc u a o cxptrfi/" 2 - 1 ). It was obtained by the trichromatic 
image (figure 0), which was restored with a 3 X 1.5 mas beam. 
The fit was applied at pixels of 0.1 X 0.1 mas with brightness 
thresholds of 5cr at 8.4 and 15.4 GHz. Pixels below 5cr at 
2.3 GHz (and above 5cr at 8.4 and 15.4 GHz) are valid for 
the fit; we put the 5cr value as an upper limit. We clipped 
out pixels where the FFA fit did not result in a positive Tf. 
(b): The opacity profile along the jet in P.A. = 65° . It peaks at 
Tf ~ 300 in the nucleus. It tails over ~ 0.7 pc along the western 
(receding) jet, and rapidly falls at 0.1 pc in the eastern (ap- 
proaching) jet. Out of this region, an extended profile of an 
ambient FFA absorber can be seen. The green curve shows the 
best-fit model of the FFA opacity by the ambient plasma with 
the isothermal King model distribution described as equation 
(12). (c): A schematic diagram of NGC 1052. The nucleus 
and the double-sided jets are visible at 15.4 GHz, with the 
viewing angle of 50° . The ambient plasma (green) distributes 
around the nucleus supported by its thermal pressure, which 
contributes an extended FFA components. The FFA torus 
is presumed to lie perpendicular to the double-sided jets. It 
covers ~ 0.7 pc of the receding jet and ~ 0.1 pc of the jet, to 
generate the asymmetry seen in the opacity distribution. 
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T b = (2.32 ±0.04) x 10 9 K and (0.97 ±0.02) x 10 9 K, re- 
spectively, while component B is unresolved and shows the 
highest Tb > 1-48 x 10 10 K. Thus, we identify component 
B as the nucleus. 

3.2. Viewing Angle of the Jet Axis 

To clarify the geometry of the jet, we 
analyze the Doppler beaming effect (O00 
[cite]cite.Blandfordl997Blandford, K6nigl(1997); |p0 
[cite] cite. Pcarsonl987Pcarson, Zensus(1987)) . Although 
both FFA and the Doppler beaming result in an apparent 
asymmetry of brightness between approaching and re- 
ceding jets, the former is negligible at frequencies beyond 
the peak. If we assume the intrinsic symmetry of the 
double-sided jets, the brightness ratio, R, is related to 
the intrinsic jet velocity, (3 (in the unit of the speed of 
light), and the viewing angle, 9, by 



Tr 



1 + (3cos9 
1-/3 cos 9 



(3) 



Here, T b and T b arc the brightness temperature of the 
approaching and receding jets, respectively. Substituting 
T+ = (2.32 ±0.04) x 10 9 K and = (0.97 ±0.02) x 10 9 
K of components A and C, respectively, we have R = 
2.4 ±0.1 at 15 GHz, and /3cos0 = 0.16 ±0.01 for a = 
—0.65. The apparent proper motion of the jet gives 
another constraint. The relative proper motion of 1.3 
mas yr _1 between the approaching and receding jets (^00 
[cite]cite.Kellermannl999Kellermann et al.(1999)) corre- 
sponds to /3 app = + f3~ — 0.40 . Here, the positive and 



negative signs are for the approaching and receding jets, 
respectively. Using the formula of the apparent velocity 
(000 [cite]cite.Blandfordl997Blandford, K6nigl(1997)), 



ft 



app 



l + /3cos0' 



we have /3sin(? = /3 app (l — /3 2 cos 2 9)/2. Finally, we ob- 
tain (3 = 0.25 and 9 = 50°. Although this is a rough esti- 
mation, it is obvious that the eastern and western jets 
are approaching and receding, respectively. A proper- 
motion analysis based on the VLBA monitoring program 
at 15 GHz (ij00 [cite]cite.Kellermannl998Kcllermann et 
al.(1998); [DO [cite]cite.Kellermannl999Kellermann et 
al.(1999)) will confirm these results and evaluate the ac- 
curacy. 

3.3. The Absorption Mechanisms 

We discuss the origin of the low-frequency cutoff 
in the spectra: SSA or FFA. In the SSA process, 
the opacity is related to the magnetic field, B, and 
brightness temperature, Tb, at the peak frequency. If 
we assume equipartition between the magnetic field 
and the electrons, the peak frequency, u m , in GHz is 
given by [t]e]cite.Kellermannl981Kellermann, Pauliny- 
Toth( 1 98 1 ) ( [cite] cite . Kellermannl 98 1 Kellermann, 

Pauliny-Toth(1981)), 
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Fig. 5. Tolerance of the FFA opacity profile along the jet 
against the position errors in registration. The thick line is 
the same as in figure which shows the opacity profile with 
nominal registration. The thin lines represents the results of 
four test cases, with shifts at 15.4 and 2.3 GHz, (<5§ 15 - 4 ,5§ 2 - 3 ) 
of (+0.2, +0.2), (+0.2, -0.2), (-0.2, +0.2), (-0.2, -0.2) mas. 
These tentative profiles align with that of the nominal regis- 
tration, representing the peak towards the nucleus, the tail to 
the west receding jet, and the break at 1 mas to the east ap- 
proaching jet. A spike of Tf = 78 at —2.7 mas from the nucleus 
can be seen only in the case of (<5£ 15 ' 4 , 5£ 2 ' 3 ) = (-0.2, +0.2) 
mas, which reflects on the tolerance of the registration when 
opposite registration errors for 2.3 and 15.4 GHz. 

Here, the magnetic field, B, is in G, the peak flux density, 
S m , is in Jy, the component size, <f>, is in mas. Equation 
(||) is replaced by 

S = 4.57x 10 19 zy m T- 2 (l + z), (6) 

using the brightness temperature of the Gaussian compo- 
nent, Tb, in K, 

T b = 1.224 x lO 12 SV~ 2 0- 2 (l + z). (7) 

For the approaching jet component, A, the peak flux den- 
sity, S m = 2.0 Jy, and the size of <f> = \J </>maj 0min = 2 mas 
at the peak frequency, v m = 4.9 GHz, give T b = 2.6 x 10 10 
K. In the co-moving frame, the peak frequency is modified 
as = 4.3 GHz by the Doppler correction, 

1-/3 cos 8 



(8) 



The brightness temperature in the co- moving frame, T b , 
is also corrected by 

/ \ 2-ao 



and T b * = 1.5 x 10 10 K is obtained. We thus have B = 0.9 
G. 

For the receding jet component, C, we have v m = 7.6 
GHz, S m = 0.25 Jy, = 1.0 mas, and T b = 5.3 x 10 9 K. 
These values should be corrected as 
l+/3cos(9 



^~8S 1 /55^V- 4/5 (l + ^) 1/5 - 



(10) 



(5) 
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Thus, we have v* m = 9.2 GHz and T* = 8.6 x 10 9 K; there- 
fore B = 5.7 G. Such strong magnetic fields arc unrealistic 
and impossible to sustain synchrotron emission for years 
against synchrotron loss. If we assume that components A 
and C have been emitting for 20 yr (= separation / appar- 
ent speed) without re-acceleration, B < 55 mG is required 
to sustain synchrotron emission with a break frequency 
of 1st > 15.4 GHz. In this case, v m < 260 MHz. Thus, 
the synchrotron emission is likely to be optically thin, in 
terms of SSA, towards the jets throughout our observing 
frequencies. 

Since SSA is an internal process in the jet, a difference in 
the peak frequencies at the rest frame requires an intrinsic 
asymmetry of the magnetic fields between the approaching 
and receding jets. On the contrary, FFA naturally gener- 
ates a higher peak frequency towards the receding jet than 
towards the approaching one. In the case of a spherically 
symmetric distribution of the external absorber, the path 
length along the line of sight will be longer towards the 
receding jet than towards the approaching one (see figure 
. Furthermore, if the absorber forms a disk or a torus 
perpendicular to the jet, it covers the receding jet, but 
does not cover the approaching jet. Such a geometry is 
discussed in the next section. 

Although Tb and v m are uncertain for component B, 
the magnetic field can be roughly estimated. Since the 
total flux density of 1.02 ±0.12 Jy at 31.4 GHz (§00 
[cite]cite.Geldzahlerl981Geldzahler, Witzel(1981)) is less 
than a natural extension of the rising spectrum of com- 
ponent B, the peak frequency, f m , should be 15.4 GHz 
< v m < 31.4 GHz. The intrinsic brightness temperature 
should be less than the Compton limit, Tb < 10 12 K. 
Therefore, 7 x 10 -4 G < B < 7 G. Although the nucleus 
component, B, could be optically thick by SSA, its spec- 
tral index of oig 5 ^ 4 = 3.1 ± 0.1 requires a significant con- 
tribution of FFA. To keep simplicity, we later analyze the 
spectra assuming that the synchrotron emission is opti- 
cally thin in terms of SSA, rather than a mixture of SSA 
and FFA. Our FFA analysis would be modified if the SSA 
opacity were revealed by future observations. 

In some GPS sources, the SSA model out- 
puts an acceptable strength of the mag- 
netic fields. §c]cite.Mutell985Mutel et 
al.(1995) ([cite]cite.Mutell985Mutel et al.(1995)) 
calculated magnetic fields of ~ 10~ 4 — 10~ 3 G 
from SSA for three sources: 1518+047, CTD 93, 
and 2050+364. [t]e]cite.Readheadl996Readhead et 
al.(1996) ([citc]citc.Readhcadl996Readhead et al.(1996)) 
estimated magnetic fields of 3 x 10~ 2 /i 2 / 7 G and 
4 x 10~ 2 /i 2 / 7 for the bright jet components, B\ and B2, 
of a GPS source 2352+495, based on the brightness 
temperatures of T b = 2 x 10 10 K and 2 x 10 11 K, and a 
peak frequency of v m ~ 1.66 GHz, with a Hubble constant, 
H , of 100ft km s" 1 Mpc" 1 . §3]cite.Shafferl999Shaffer 
et al.(1999) ([cite]cite.Shafferl999Shaffer et al.(1999)) 
derived B = 10~ 4 G for the lobes of a GPS source, CTD 




Offset along the jet [pc] 

Fig. 6. FFA opacities of NGC 1052 out of the torus. The 
horizontal axis indicate the offset along the jet in P.A.= 65° 
in pc, where 1 mas corresponds to 0.1 pc. The nega- 
tive and positive offsets stand for the eastern (approaching) 
and western (receding) jets, respectively. The vertical axis 
shows the FFA coefficient, Tf, derived from the FFA model 
fit. The opacities in the torus, —1 < x < 4 mas (or -0.1 
to 0.4 pc), were removed from the plot. The solid grey 
line shows the opacity calculated from the isothermal King 
model, which is calculated in the Appendix and shown in 
figure pi Here, we take the viewing angle 9 = 50°. The 
best fit gives n = (1.00 ±0.03) X 10 4 (T e /10 4 K) 3 / 4 cm" 3 and 
r c = 0.84 ±0.02 pc. 

93, based on T h = 1.3 x 10 11 K and 0.7 x 10 11 K for the 
northern and southern lobes, respectively, and u m ~ 1 
GHz. A higher Tb and a lower u m than those of NGC 
1052 result in the much lower B required from the SSA 
mod el for these sources. 

te]cite.Shafferl999Shaffer et 
al.(1999) ([cite]cite.Shafferl999Shaffer et al.(1999)) 
claimed that the observed dependence of the apparent 
component separation of CTD 93 on wavelength is 
additional evidence for SSA. NGC 1052 also shows such 
a frequency dependence of the apparent component 
separation (see table 1); nevertheless, we mention that 
the central condensation of the FFA opacities can explain 
it. 

3. 4- Free-Free Absorption Opacity Distribution 

By applying a least-squares fit for the FFA spectrum 
at every image pixel, we obtain the distribution of the 
FFA opacity, Tf , shown in figure [|a. The fit is applied for 
restored images at all three frequencies with a 3.0 x 1.5 
mas beam in P.A.= 0°. The opacity map is clipped by the 
thresholds for the brightness of 5a at 8.4 and 15.4 GHz. 
Pixels below 5<r at 2.3 GHz (and above 5<t at 8.4 and 
15.4 GHz) are valid for the fit; we put the 5<r value as an 
upper limit. We reject pixels where the FFA fit does not 
result in a positive Tf. The intrinsic spectral index, ao, is 
fixed to be —0.65 at every pixel, since three frequencies 
are insufficient to allow any freedom in ao- It should be 
mentioned that the FFA opacity would be overestimated 
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if ao were larger in the nucleus. 

Figure [|b shows a slice of the distribution of the FFA 
opacities along the jet in P.A.= 65°. The error bars in- 
dicate formal errors of the least squares fit for the FFA 
spectrum. The errors of registration between images at 
different frequencies, derived in subsection 3.1, can con- 
tribute to additional uncertainties of the FFA opacities. 
To evaluate the uncertainties responding to the registra- 
tion errors, we tentatively shifted images at 2.3 and 15.4 
GHz relative to the 8.4-GHz image by 0.2 mas, performed 
a least-squares FFA fit, and then plotted the opacity pro- 
files along the jet. The results of four test cases, with the 
shifts at 15.4 and 2.3 GHz, (<5£ 15 - 4 ,<5£ 2 - 3 ) of (+0.2, +0.2), 
(+0.2, -0.2), (-0.2, +0.2), (-0.2, -0.2) mas, are plotted 
in figure 0. The profiles align with that of the nominal reg- 
istration, except for a spike of Tf = 78 at —2.7 mas from 
the nucleus in the case of (S£_ 15A ,S£_ 2 - 3 ) = (-0.2, +0.2) mas. 
This spurious spike reflects on the tolerance of the regis- 
tration, where the brightness distribution gaps between 
components A and B. However, we have a small joint 
probability to meet such a bad case with opposite position 
errors for 2.3 and 15.4 GHz. We also attempted position 
shifts in the north-south direction, and confirmed that the 
opacity profile is insensitive in this direction, probably due 
to elongation of the restoring beam. 

As shown in equation (0), Tf can be used to investigate 
the cold dense plasma. The opacity peaks at Tf ~ 300 
towards the nucleus and decreases to both downstream 
regions. In the central sub-pc from the nucleus, we found 
a significant asymmetric distribution along the jet. The 
opacity rapidly falls at 0.1 pc in the eastern jet, while it 
smoothly decreases over 0.7 pc in the western jet. This 
feature can be seen in the trichromatic image (figure §), 
in which blue color indicates an opaque area at low fre- 
quency due to FFA. As stated before, we find that the 
eastern and western jets are approaching and receding, 
respectively, in the viewing angle of 9 ~ 50°. The dense 
absorber does not cover the approaching jet, but does 
cover the receding jet. Thes implies that the absorber 
forms a disk or torus perpendicular to the jet, as illus- 
trated in figure §c. In this configuration, the receding jet 
is covered with the torus, but the approaching jet isn't. A 
rapid change of the opacity at 0.1 pc in the approaching 
jet suggests that the torus is geometrically thick with a 
thickness H > rcotd, where r is the radius of the torus. 
Since the opacity tails along the receding jet up to ~ 0.7 
pc from the nucleus, we consider that the radius of the 
torus is ~ 0.5 pc. The maximum opacity, Tf = 300, for a 
path length of ~ 0.7 pc gives a constraint for n c and T e 
by n 2 ! 7 " 1 - 5 = 930. The ionization condition of T c > 10 4 
K constrains n c > 3.1 x 10 4 cm~ 3 , and thus the column 
density, n c L ~ 0.7 x 10 23 (T C /10 4 K) 3 / 4 cm- 2 . This is con- 
sistent with an atomic column density of ~ 10 23 cm -2 
derived from the ROSAT and ASCA X-ray observations 
(§00 [cite]cite.Guainazzil999Guainazzi, Antonelli(1999)). 
Such a high column density is found via H I absorption to- 
wards some Seyfert and radio galaxies. Although the Hi 
column densities hold an ambiguity of the spin tempera- 
ture, T spin , those of (3.9 + 0.5) x 10 21 (T spin /100K) cm" 2 



towards NGC 4151 (§00 [cite]cite.Mundelll995Mundell et 
al.(1995)), (2.54 + 0.44) x 10 19 T spin cm" 2 towards Cygnus 
A (000 [cite]cite.Conwayl995Conway, Blanco(1995)), 
14.0 x 10 21 (r spin /100 K) cm" 2 towards the core of 
Hydra A (§D0 [cite]cite.Taylorl996Taylor(1996)), 27.8 x 
10 23 (T spin /8000K) cm" 2 in the broad velocity-width 
component towards the core of PKS 2322—123 (§00 
[cite]cite.Taylorl999Taylor et al.(1999)), and (19.3 ± 
3.8) x 10 21 (T spin /100K) cm" 2 towards NGC 3079 (§00 
[cite]cite.Sawada2000Sawada-Satoh et al.(2000)) are com- 
parable with the electron column densities probed by FFA. 
A comparison of distribution of FFA with that of H I ab- 
sorption could illustrate the ionizing surface of the torus. 

3.5. Ambient Free-Free Absorbing Plasma 

Besides the asymmetric condensation in the central 0.1 
pc, the opacity widely extends in the outer region, x < 
— 1.5 mas, in the approaching jet and x > 10 mas in the 
receding one. Here, x is the offset from the nucleus along 
the jet in P.A.= 65°. The opacity profile in this ambient 
absorber can be explained by the isothermal King model 
(§00 [cite]cite.Kingl972King(1972)), 



n 



l+l- 



(12) 



where n is the electron density at the centroid, r the 
radius from the centroid and r c the core radius. Making a 
least-squares fit for the opacities in x < — 1 and x > 4 mas, 
we estimated n = (1.00 + 0.03) x 10 4 (T C /10 4 K) 3 / 4 cm- 3 
and r c = 0.84 + 0.02 pc. The best- fit curve is shown as the 
green line in figures and §. See the Appendix for the 
analytic profile of the opacities towards the slant jets. 

It is remarkable that the torus-like plasma distribu- 
tion coexists with the spherically symmetric distribu- 
tion. The former should be supported by the centrifu- 
gal force by rotation around the nucleus, while the lat- 
ter is supported by the thermal pressure. We consider 
that our results witness the transition stage from the ra- 
dial accretion to rotational accretion. We note that water 
masers (§00 [citc]citc.Claussenl998Claussen et al.(1998)) 
distribute between components B and C, where the FFA 
torus covers the foreground of the receding jet. Our results 
support the idea that the masers are amplified radio con- 
tinuum emission from the receding jet by the foreground 
disk. The velocity gradient along the jet axis, faster re- 
ceding velocity at a position closer to the nucleus, may 
represent a dynamical infall in the torus. 

4. Conclusion 

We have performed trichromatic VLBA observations 
towards the nucleus of NGC 1052. Based on the peak 
frequency and brightness temperature, we identified the 
nucleus component in this object. The kinematics and 
orientation of the double-sided jets were derived from the 
Doppler beaming effect. We found that the spectra of 
the core and double-sided jets can be better fitted with 
the free-free absorption process by the cold dense plasma. 



No. 2] 



Plasma Torus around NGC 1052's Nucleus 



9 



The distribution of the opacity is biased towards the re- 
ceding jet, suggesting a torus-like condensation perpen- 
dicular to the jet. Out of the torus, we also found an 
ambient plasma distribution which follows the spherical 
isothermal King model. Our results indicate the first look 
at a dynamical transition from thermally-supported into 
centrifugally-supported in the subparsec-scale region of 
AGN. 

We thank M. Miyoshi for helpful discussions. The 
VLBA is operated by the National Radio Astronomy 
Observatory, which is a facility of the National Science 
Foundation (NSF) operated under a cooperative agree- 
ment by Associated Universities, Inc. The Green Bank 
Interferometer is a facility of the NSF operated by the 
NRAO in support of NASA High Energy Astrophysics 
programs. This research has made use of data from the 
UMRAO, which is supported by funds from the University 
of Michigan. 



Appendix. FFA Opacities along the 
Sphere with the Isothermal King Model 



Plasma 



Here, we discuss the optical depth of free-free 
absorption in the plasma sphere, whose electron 
density follows the isothermal King model (QOO 
[cite]citc.Kingl972King(1972)). The model is derived 
from the force balance between the pressure gradient and 
the self-gravity of the gas, assuming spherical symmetry 
and thermal equilibrium. The model gives the radial pro- 
file of the particle density, p, as a function of the radius, 
r, by 



receding jet 




Fig. 7. FFA optical depth towards the twin jets along the 
plasma with the density profile of the isothermal King model. 
Let the xi and x axes be parallel and perpendicular to the 
line of sight, respectively, and let the nucleus be the origin. 
The jet, as a synchrotron emitter, is inclined from the line of 
sight by the viewing angle 9. The opacity integration, as a 
function of x, is taken from the jet at zq = — ta ^g to infinity. 



p(r) = pa 



l+l- 



(Al) 



Here, po is the density at the center and r c is the core ra- 
dius. If we assume the isothermal condition, the relation 
between the electron density, n c , the electron tempera- 
ture, T c , and the FFA opacity, tq, will be 



T = 0.46 



T r 2 nddL 



LOS 



0.46T C 



0.46T; 



nidL 



LOS 



'cO 
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1 + 


(0} 


/LOS 







dL. 



(A2) 



and 



Here, n e rj is the electron density at the center in cm" 
L is the path length along the line of sight in pc. 

Let us consider the configuration as shown in fi gure 
[?]. The integration of the optical depth in equation ( [A 2] ) 
should be taken from the jet at £o = — tl^T? to infinity. The 
radius from the nucleus, r, is describe d in (x,£) coordinate 
+ £ 2 . Therefore, equation (A2) will be 



r (x) =0.46r; 



; 2 6 
n cO r c 



(A3) 



Here, we use the mathematical formula 



Isothermal King model 




-0.5 0.5 1 

Offset along the jet [in unit of the core radius] 

Fig. 8. Profiles of the FFA optical depth towards the jets 

with different viewing angles. The horizontal axis is the offset 

along the jet in unit of the core radius, r c . The negative offset 

stands for the approaching jet side. The vertical axis is the 
_ 3 

normalized opacity, when T e 2 n^ r c is unity. 



10 



S. Kameno et al. 



[Vol. 53, 



(x 2 + C)~ 3 dx 



4C*(x 2 + 
3 1 

sc 2 ^ 
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arctan 
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8C 2 (x 2 - 
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Figure || shows opacity profiles along the jets wit h dif- 
ferent viewing angles, calculated by equation (|A6|). As [Readhead et al.(1996)] Readhead, A. C. S., Taylor, G. B., Xu 



W., Pearson, T. J., Wilkinson, P. N., & Polatidis, A. G. 
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enough from the center, the opacity is small at any given [Shepherd(1997)] Shepherd, M. C. 1997, in Astronomical Data 



easily supposed, the slant jet axis causes an asymmetric 
profile of the opacity. Since the electron density is larger 
towards the center, the opacity is larger within the core 
radius than outside. When the radiation source is far 



viewing angle 
References 

[Asada et al.(2000)] Asada, K., Kameno, S., Inoue, M., 
Shen, Z.-Q., Horiuchi, S., & Gabuzda, D. C. 2000, 
in Astrophysical phenomena revealed by Space VLBI, 
ed. H. Hirabayashi, P. G. Edwards, & D. W. Murphy, 
(Sagamihara, ISAS), 51 

[Bicknell et al.(1997)] Bicknell, G. V., Dopita, M. A., & O'Dea, 
C. P. 1997, ApJ, 485, 112 

[Blandford, K6nigl(1997)] Blandford, R.D., & Konigl, A. 1979, 
ApJ, 232, 34 

[Claussen et al.(1998)] Claussen, M. J., Diamond, P. J., Braatz, 
J. A., Wilson, A. S., & Henkel, C. 1998, ApJ, 500, L129 

[Conway, Blanco(1995)] Conway, J. E., & Blanco, P. R. 1995, 
ApJ, 449, L131 

[de Vries et al.(1997)] de Vries, W. H., Barthel, P. D., & O'Dea, 

C. P. 1997, A&A, 321, 105 
[Geldzahler, Witzel(1981)] Geldzahler, B. J., & Witzel, A. 

1981, AJ, 86, 1306 
[Guainazzi, Antonelli(1999)] Guainazzi, M., & Antonelli, L. A. 

1999, MNRAS, 304, L15 
[Ho et al.(1997)] Ho, L. C, Filippenko, A. V., Sargent, W. L. 

W., & Peng, C. Y. 1997, ApJS, 112, 391 



Analysis Software and Systems VI, ed. G. Hunt, & H. E. 
Payne, (San Francisco, A. S. P. Conference Series), 77 
[Taylor(1996)] Taylor, G. B. 1996, ApJ, 470, 394 
[Taylor et al.(1999)] Taylor, G. B., O'Dea, CP, Peck, A. B., & 

Koekemoer, A. M. 1999, ApJ, 512, L27 
[Vermeulen et al.(1994)] Vermeulen, R. C, Readhead, A. C. S., 

& Backer, D. C. 1994, ApJ, 430, L41 
[Walker et al.(2000)] Walker, R. C, Dhawan, V., Romney, J. 
D., Kellermann, K. I., & Vermeulen, R. C. 2000, ApJ, 530, 
233 

[Walker et al.(1994)] Walker, R. C, Romney, J. D., & Benson, 
J. M. 1994, ApJ, 430, L45 



